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Abstract 

We define the notion of the Fourier transform for the rook monoid (also 
called the symmetric inverse semigroup) and provide two efficient divide- 
and-conquer algorithms (fast Fourier transforms, or FFTs) for computing 
it. This paper marks the first extension of group FFTs to non-group 
semigroups. 

1 Introduction 

The rook monoid also called the symmetric inverse semigroup, is the set 
of all partial permutations of {1, . . . ,n} under function composition, i.e., it is 
the semigroup of all bijections between all subsets of {1, . . . , n} under function 
composition, with the usual rule for composing partial functions. That is, gof is 
defined for precisely the elements x for which x € dom(/) and f{x) G dom(g). 
It is called the rook monoid because it is isomorphic to the semigroup of all 
n X n matrices with the property that at most one entry in each row is 1 and at 
most one entry in each column is 1 (the rest being 0) under multiplication. Such 
matrices (called rook matrices) correspond to the set of all possible placements of 
non-attacking rooks on an n x n chessboard. For example, consider the element 
(7 G i?4 defined by 

cr(2) 1, cr(4)=4. 

Then, viewed as a partial permutation, a is 



12 3 4 

-1-4 



where the dash indicates that the above entry is not mapped to anything. As a 
rook matrix, we have 

"0100 


1 
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In this paper we define the notion of the Fourier transform for a C-valued 
function on the rook monoid, and we address the problem of computing it 
efficiently. We present two distinct fast Fourier transform (FFT) algorithms, 
both of which make use of tools generalized from group FFT theory [25] . We 
assume no specialized knowledge. This paper is meant to be readable by group 
FFT theorists with an interest in semigroups as well as semigroup theorists with 
an interest in FFTs. 

The history of the FFT is an interesting one, beginning, like many subjects 
in mathematics, with Gauss, who developed the now-classical FFT to efficiently 
interpolate the orbits of asteroids. This FFT was rediscovered by Cooley and 
Tukey [5] in 1965 to analyze seismic data. A different FFT was discovered in 
1937 by Yates [33j, for the analysis of data on factorial designs. The theory 
of FFTs for Abelian groups reconciles these two FFTs. The first is an FFT 
on Z/nZ, and the second is an FFT on (Z/2Z)'^. From this generalization also 
grew a theory of FFTs for non- Abelian groups. For example, we now have FFTs 
on supersolvable groups, FFTs on the symmetric group, FFTs on compact Lie 
groups, and FFTs on finite groups of Lie type. For a survey of these results, see 

m- 

We have both theoretical and practical motivations for creating FFT algo- 
rithms for the rook monoid. From a theoretical perspective, these FFTs are 
significant because they are the first examples of FFTs on a non-group semi- 
group. The rook monoid plays the same role for finite inverse semigroups that 
the symmetric group Sn does for finite groups in terms of Cayley's theorem: 
every finite inverse semigroup is isomorphic to a sub-semigroup of some rook 
monoid (see, e.g. [17], p. 36-37), so this makes the rook monoid a natural place 
to begin extending the theory of FFTs to semigroups. Also, just as groups cap- 
ture symmetry, inverse semigroups capture partial symmetry. See [16] for more 
on this idea. 

From a practical perspective, partially ranked data occurs naturally, and the 
rook monoid provides a nice computational framework for analyzing such data. 
For example, consider voting data in which some voters do not fully rank all 
candidates. In practice, this may well be the case. For example, in the 1980 
American Psychological Association (APA) presidential election, only a third of 
the voters cast fully ranked ballots [8]. Say f{a) is the number of voters who 
submitted a ballot of type a. For the a above, f{a) would represent the number 
of voters who ranked candidate 2 in the first position, candidate 4 in the fourth 
position, and didn't bother ranking candidates 1 and 3. The values {/(c)} make 
up the dataset. There have been previous attempts at spectral analysis for such 
datasets. For example, in [8], Diaconis used spectral analysis on symmetric 
groups to analyze the APA election mentioned above. We propose viewing / as 
an element of the semigroup algebra Ci?„, so that spectral analysis of / involves 
decomposing Ci?„ into its isotypic components by means of an FFT. Things are 
more complicated than in the group case because Ci?„ has two natural bases. 
This means that there are two natural ways to view functions on i?„ as elements 
of CRn, and therefore two different notions for convolution of such functions. 
We will address both in this paper. 
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Partially ranked voting data is also a type of "not missing at random" 
(NMAR) data. Missing data in a sample is said to be NMAR if it is believed 
to be missing, at least in part, because of its unobserved value. In the voting 
example, any ballot that is not fully ranked is considered NMAR, since the 
voter intentionally left something on the ballot blank. Consumer survey data 
frequently contains NMAR data. NMAR data contrasts with randomly missing 
data, which, in the voting example, can be caused by accidentally losing some of 
the ballots. To date, there doesn't seem to be a standard method for analyzing 
NMAR data. However, the concept of spectral analysis for the rook monoid 
immediately applies to any NMAR dataset that can be considered a collection 
of partial rankings. Since spectral analysis is a model-independent approach, 
spectral analysis for the rook monoid may go a long way towards creating a 
standard method for dealing with NMAR data. 

This paper is organized as follows. Sections 12.11 and [ 2.21 contain basic defini- 
tions and background material about semigroups and the rook monoid. Section 
12.31 contains a brief introduction to semigroup representation theory. Sections 
12.41 and 12.51 are devoted, respectively, to two important tools in FFT theory: 
chain-adapted matrix representations and Schur's Lemma, adapted to semi- 
groups. 

Section [3] is devoted to basis considerations for semigroup algebras. Section 
13.21 defines both natural bases for Ci?„, and Section [231 defines the notion of a 
Fourier basis. Section [Jj defines the Fourier transform. 

Section [5] shows how the semigroup algebra Ci?„ decomposes into a product 
of matrix algebras over group algebras. This has far-reaching implications, 
some of which are explored in Section [6] Section 16.11 describes explicit matrix 
representations for Section [6.21 contains an explicit Fourier basis description 
for CRn, and Section contains the Fourier inversion theorem for Ci?„. 

In Section [71 we explain some of the issues that arise when choosing how to 
associate functions on i?„ to elements of Ci?„. As mentioned previously, this is 
only an issue because Ci?„ has two natural bases. 

Section [51 contains FFT-specific considerations. Computational complexity 
for FFTs is defined in Section [5TT1 Our algorithms for the FFT on i?„ use ideas 
from group FFT theory, which we cover in Sections 18.21 and 18.31 

Section [9l contains our better FFT algorithm for i?„ . It relies heavily on the 
decomposition described in Section \5\ and the results in Sections 16.11 and 18.31 

Section [TOl contains our other FFT algorithm for i?„ . The ideas involved 
are mainly from Sections 12. 4( 12.51 and 18.21 While it is not as efficient as the 
algorithm given in Section [SI it is constructed in an entirely different way. We 
present this algorithm because the ideas involved may be useful for constructing 
FFTs on other semigroups. 
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2 Preliminaries 



2.1 Definitions for general semigroups 

Definition (semigroup). A semigroup is a nonempty set S together with an 
associative binary operation {which we write multiplicatively) . 

Definition (inverse semigroup). An inverse semigroup is a semigroup S such 
that, for each x € S, there is a unique y d S such that 

xyx — X and yxy — y. 

In this case, we write y = x^^ . 

We remark that the condition that y be unique is necessary for this definition. 
An element x G 5* is said to be regular or Von-Neumann regular if there is at 
least one y (z S satisfying xyx = x and yxy = y, and S is said to be regular 
if every element of S is regular. Consider the full transformation semigroup 
X on the set {1, 2, . . . , n}; that is, all maps from {1, 2, . . . , n} to itself under 
composition. It is easy to see that X is regular, and that (for n > 2) there exist 
elements x E X for which there are multiple elements y € X satisfying xyx — x 
and yxy = y. X \s therefore not inverse. An equivalent characterization of 
inverse semigroups (see, e.g. [17]) is as follows. 

Definition (inverse semigroup). An inverse semigroup is a semigroup S which 
is regular and for which all idempotents of S commute. 

Definition (semigroup algebra). Let S be a finite semigroup. The semigroup 
algebra CS is the formal C-span of the symbols {s}ses- Multiplication in CS, 
denoted by *, is given by convolution {i.e., the linear extension of the semigroup 
operation via the distributive law): Suppose /, g £ CS, with 

res teS 

Then 

/*5 = E/w^E5Wt-E E f(^)9(t>- (1) 

res teS seS r,teS:rt=s 

If S' is a group, then convolution may be written in the familiar way: 

/ * 5 = E E f(^)9{r~'^s)s. 
seSreS 

Note that convolution here means convolution in the semigroup algebra. We 
hold off on defining a notion of convolution for functions on S until Section [71 
We will use the semigroup algebra convolution to define function convolution 
once we have defined how to associate functions to elements of the semigroup 
algebra. It turns out that, unlike in the case of groups, there is in general more 
than one natural way to associate functions on S to elements of CS because, in 
general, CS has more than one natural basis. 
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2.2 Facts about 

Definition (rank). Given an element a G Rn, the rank of a, denoted rk((T), is 
defined to be rk((T) — |dom((T)| = |ran(cr)|. 

There are two main types of elements of i?„: the elements of rank n (i.e., 
the permutations) and the elements of rank less than n. In his analysis of the 
representation theory of i?„, Munn [Mj introduced what he called cycle-link 
notation for the elements of Rn ■ As an example, consider the element a ^ R4 

12 3 4 

3-24 

In cycle-link notation, a cycle (ai, 02, . . . , a^) means that 

ai 1-^ 02, 02 1-^ 03, ... , afe„i i-> Ok, and Ofc 1-^ ai, 

and a link [61, 62, ... , 6^] means that 

&i ^ 621 ^2 ^ . . . 7 &fc-i ^fci and 6fc goes nowhere. 

The element a, expressed in cycle-link notation, would be [1, 3, 2] (4). Note that 
the cycle-link representation of a permutation always consists of cycles only, 
and elements of rank less than n always contain links when written in cycle-link 
notation. 



Theorem 2.1. 



n / \ 2 

k=0 ^ ' 



Proof. For any particular rank k, there are (^) choices for the domain and (^) 
choices for the range of an element of i?„ , and for any particular choice of domain 
and range, there are k\ ways of mapping the domain to the range. □ 

We also have the recursive formula 

Theorem 2.2. For n > 3, 

\Rn\ = 2n\Rn-i\ - (n - l)^\Rn-2\- (2) 
Proof See Section [THSl □ 
Theorem 2.3. i?„ is an inverse semigroup. 
Proof For an element a £ Rn, define 7 £ i?„ by 

dom(7) = ran(o'), and, for x G dom(7), 

7(x) = (T^^{x). 

It is easy to see that 7(77 — 7 and a^a — a, and that 7 is the only element of 
Rn satisfying both equations. □ 
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2.3 Semigroup Representations 

Here we provide a brief introduction to the representation theory of semigroups. 
For a more complete introduction to this topic, see [53] and [5T| . 

In this paper, we adopt the convention that all maps act on the left of sets, 
and we only consider left modules. We also only concern ourselves with finite 
semigroups. For the rest of the paper, every semigroup discussed is understood 
to be finite. Let 5 be a semigroup. If S does not have identity, then let 
denote S with an identity element added. If S does have identity, then let 
5*1 = S. 

Definition (semigroup representation). A representation of S is a finite di- 
mensional C-vector space V together with a homomorphism p from S into the 
semigroup End(V^) under composition. 

Remark: If 5 is a group, this definition docs not agree with the usual group 
definition, which requires p{e) — 1 G End(V^). Such representations are called 
unital. However, there is only a trivial difference between a representation and 
a unital representation of S*^: a representation is called null if p{s) — for all 
s S, and every representation of is either unital, null, or a direct sum of a 
unital and a null representation ([H], Fact 1.10). 

Definition (semigroup matrix representation). A matrix representation of S 
is a representation of S with respect to a particular choice of basis for V , i.e., 
it is a homomorphism from S into the semigroup of Ami{V) x dim(V^) matrices 
with entries in C under multiplication. 

Definition (representation module). A representation module for S is a CS- 
module that is also a finite dimensional C-vector space. 

Definition (matrix representation of CS). A matrix representation of the semi- 
group algebra CS is an algebra homomorphism p from CS into the matrix algebra 
Md,{C). 

Given a representation p of S*, we may create a representation module for 
S by viewing the underlying vector space F as a C5'-module, where the action 
of 5' on y is given by p, and the action of CS' on V is given by the linear 
extension of p. In this case, we call V (as a CS-module) the representation 
module for p. Conversely, given a CS*- module V which is also a finite dimensional 
C-vector space, we may define a representation p : S ^ End(l^) by p{s) ■ v 
= s ■ V for all s e S* and v d V. In this case, we call p the representation 
associated to V. Furthermore, a matrix representation of CS defines a matrix 
representation of S by restricting to a basis, and any matrix representation of 
S extends linearly to a matrix representation of CS. Matrix representations 
of S and of CS arc therefore in one-to-one correspondence, and the notions of 
representations, matrix representations, and representation modules for S are 
equivalent. We shall move between these notions freely as needed. 

Also, we will often abuse the language and simply refer to a representation or 
matrix representation p of S, where the underlying vector space V is understood. 
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Definition (irreducible representation). A non-null representation p of S is 
said to be irreducible if the representation module for p is non-null and simple. 
In other words, with respect to no basis does p have the form 

' bi 
9 <p2 

for some representations (j)i, (j)2 and some matrix- valued function g. 

A matrix representation of CS* is irreducible if and only if the corresponding 
matrix representation of S is irreducible. 

Definition (equivalence of representations). Given two matrix representations 
Pi and p2 of S {or ofCS), pi is said to be equivalent to p2 if there exists an 
invertible matrix A such that Api(s)A^^ = P2{s) for all s G S. Equivalently, 
two representations of S are equivalent if their representation modules are iso- 
morphic (as <CS -modules) . 

Theorem 2.4 (decomposing representations). For any S such that CS is semi- 
simple, any representation of S {or of CS) is equivalent to a direct sum of 
irreducible and null representations of S {resp. CS). Furthermore, there are 
only finitely many inequivalent, irreducible representations of S {resp. CS). 

Proof. This is 25J, Theorem 1.18 and Proposition 1.19. □ 

Theorem 2.5 (Wedderburn-Artin). Let CS be semisimple, and let y be a com- 
plete set of inequivalent, irreducible matrix representations for S. For each 
p G y, p : S ^ Md^ {C); extend p linearly to p : CS Md^ (C). Then the family 
{p} defines an isomorphism of algebras 

C5^0A/,^(C). (3) 

Explicitly, 

E - ^ (e - (e ■ 

ses p&y \s6S / pi^y \s&s / 

Proof. See, for example, [7], [25], [28], or [11], p. 50, exercise 18. □ 

Corollary 2.6 (sum of the squares of the dimensions). Let S be a semigroup 
such that CS is semisimple, and let y be a complete set of inequivalent, irre- 
ducible matrix representations for S . Then 

1^1 = E< (4) 

pey 

Proof. The formula ([4]) is just the C-dimcnsionality of the algebras appearing 
in ©. □ 
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Theorem 2.7 (Munn [53], Theorem 3.1). Ci?„ is semisimple. 

More generally, we have the following extension of Maschke's theorem. 

Theorem 2.8 (Munn [24] . Theorem 4.4). If S is an inverse semigroup, then 
CS is semisimple. 

2.4 Adapted Representations 

Let Xq < Xi < . . . < Xn be a chain of semigroups whose semigroup algebras 
<CXi are semisimple. Even if a matrix representation of Xi is irreducible, its 
restriction to will typically not be. However, it will be equivalent to 

(though not necessarily equal to) a direct sum of irreducible and null matrix 
representations of Xi^i. 

Definition (adapted representations). Let yi be a set of inequivalent, irre- 
ducible matrix representations for Xi. The collection {3^i}"=o ■^'^^'^ chain- 
adapted to the chain Xq < Xi < . . . < Xn if, for every i>l and every p G yi, 
p\xi^i is equal to a direct sum of null representations and matrix representa- 
tions in yi-i. Note that this definition forces the irreducible representations 
appearing in such a restriction to be equal whenever they are equivalent. 

Induction shows that a representation from a chain-adapted set of represen- 
tations may be restricted further down the chain with the same equality results. 
In particular, if {3^i}"^o chain-adapted to Xq < Xi < . . . < Xn, then {yj}jej 
is chain adapted to <j£j Xj for any nonempty subset J C {1, . . . , n}. 

We shall often simply ask that a complete set of irreducible representations 
yn for Xn be adapted to the chain Xq < Xi < . . . < Xn- In this case, the 
choices of the are understood (and, in fact, are often completely determined) 
by the choice of 3^„ . 

Remark: Adapted representations are very important in the construction of 
FFTs (see ^7}, for example), but the requirement that a set of representations 
be adapted is in no way limiting. Given any chain Xq < Xi < . . . < Xn 
of semigroups whose semigroup algebras CXi are semisimple, a straightforward 
induction argument shows that a complete set of inequivalent, irreducible, chain- 
adapted matrix representations always exists. Explicitly describing such sets, 
however, is frequently a challenging endeavor. 

Remark: Chain-adapted representations are also sometimes referred to as semi- 
normal or Gelfand-Tsetlin representations. 

2.5 Schur's Lemma 

Schur's Lemma appears in many forms, the most widely recognized probably 
being the one below. 
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Lemma 2.9 (Schur's Lemma). Let R be a ring and let M,N be simple R- 
modules. Then every R-module homomorphism (f) : M ^ N is either or an 
isomorphism. See, e.g. Jllf , p. 30-32. 

For computational purposes, the following (also called Schur's Lemma) is a 
very useful consequence of Lemma 12.91 

Lemma 2.10 (Schur's Lemma). Suppose that A < B < C are semigroups, that 
yA, ys, and yc are complete sets of inequivalent, irreducible matrix represen- 
tations for A, B, C respectively, adapted to the chain, and that <CA, CB, and CC 
are semisimple. Let p G yc . Say p\b — ® ■ ■ ■ ® p*' {with each p' either in ys 
or, without loss of generality, null of dimension i) and p-* |a = pi © • • • ffi P'g(j) 

{with each pf either inyA or, again without loss of generality, null of dimension 
1). Let (7 E B such that a commutes with A. Then p{a) is a block matrix 

... \ 

P^{(t) ... 

■■■ ■•• : 







pM = 







V 









p''{a) J 



where, for I < j < k, the block p^ {a) is itself a block matrix, with blocks indexed 
by all ordered pairs from {p{, . . . , Pg(j)}-' 



p'{a) 



Pi 
Pi 



Pi 

Ail/ 



Pi 

^2,2 



where the ^ are scalars, L is the appropriately-sized identity matrix for the 
block in which it appears, and the block in position p^^ , is non-zero only if 
and are equivalent representations of A. 

Proof. To prove this, we need only show that p^ {a) has the block form indicated 
above. Let be the representation module for p' on B. Viewing as a €.A- 
module, we have = (B ■ ■ ■ ® ^g{j) ' ^^'^ the Vf are simple CA-modules 

where the module action on Vf is given by . Note that a basis B of has 
already been chosen by virtue of the fact that p^ is in matrix form, and likewise, 
bases Bi for the are given since pi is in matrix form. Since p^ is adapted, we 
have that Bi C B for all i, Bi Ci Bk = 9 tor i ^ k, the Bi are ordered as subsets 
of 6, and UBi = B. 

Since cr commutes with A, we have that, for all v € and x € A, 



pi {a) 



X ■ V 



fp {ax) 



p? (xcr) ■ V — X ■ fp (cr) • w. 
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i.e., p-'icr) is a CA-lincar map from to itself. Hence 

p\a) e RomcA{V\V^) = HomcA(T^/ © • • • © ^^/o"), 
According to our basis B for , we have 

/ HomcACl^/,V/) HomcAiVlV^) .. 



HomcA ( V/ , Vi ) HomcA ( Vj-' , Vi ) 



Let X G HomcA(Vj^ , VJ;') be given in matrix form with respect to the bases 



Ba, Bb- By Lemma 12.91 X is either or an isomorphism (i.e. X is either 
or and are equivalent representations of A). Suppose then that X is an 
isomorphism. The goal, then, is to show that X is a diagonal matrix. Since X 
is CA-linear, for every x £ CA, v e we have 

Xx ■ V — X ■ Xv. 

With respect to Ba, B^ this is 

Xpi{x)v = pi{x)Xv, 

and hence for every x G CA 

Xpiix) - pi{x)X. 

Since and are equivalent, and are either null or part of an adapted set of 
matrix representations, we have = P^i and thus 

Xp^aix) = p>a{x)X (5) 

for all X e CA. 

Now, either is null of dimension 1 (in which case so is p^), which means 
that X is 1-dimensional and we're done, or p^ is an irreducible representation 
of A. So, suppose p^ is irreducible. Then by Burnside's theorem (Theorem 1.14 
of [25]), we have 

p^,(CA)=Af|H„l(C), 
and therefore ^ says that X is in the center of (C), i.e. X is diagonal. □ 

Remark/Notation: Schur's Lemma says that adapted representations cause 
the matrix p(cr) to be sparse and structured. Under the same hypotheses as 
above, let M.{B,A) be the maximum multiplicity of an irreducible or dimen- 
sion l-nuU representation of A occurring in the restriction, from B to A, of an 
irreducible representation of B. Since p{a) is a block matrix of the form 








V 











p\a) j 
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and the block [a) is of the form 



p>{a)- 



Pi 



Pi 


Pi . 






Ai^2/ • 


■ KgU)^ 


Ail/ 


Ai2 


■ ^igU)^ 



'^aU),gU)'^ J 



where ^ ^ impUes and pi are equivalent representations of A, we see that 
the matrix p{a) contains at most A4{B, A) non-zero entries per row and column. 
The computational implication of Schur's Lemma, then, is that for an arbitrary 
dp X dp matrix H , it requires no more than A4{B, A)d^ complex multiplications 
and additions to perform each of the matrix multiplications p{(j)H and Hp{a), 
as opposed to the upper bound of dp multiplications and additions that would 
be necessary if p{a) were arbitrary. 



3 Basis considerations for CS 

3.1 The poset structure of an inverse semigroup 

Definition (poset structure of S). Let S be a finite inverse semigroup. For 
s,t G S, define 

s < t <;=^ s ~ et for some idempotent e d S 
s — tf for some idempotent / G S". 

For Rn, the idempotents are the restrictions of the identity map. If S" is a 
group, then its poset structure is trivial. 

If P is a finite poset, then the zeta function C of P is given by 



C:PxP^{0,l} 

C{x,y) = 



1 if X < y 
otherwise. 



Given a poset P, one may define an incidence algebra for P over any ring 
with identity. The element C is invertible in the incidence algebra, and its inverse 
is called the Mobius function p,. There is a general theory of Mobius inversion 
for incidence algebras. We will not go into the details here. Rather, we will 
record only the results that we need. Details may be found in [31] . 

The Mobius function for i?„ (over C) is well known [30], [31]. It is 

^(x,2/) = (-l)*(-)-'-'^(^). 



11 



3.2 Natural bases for CS 

Let S be an inverse semigroup. There are two natural bases for CS (if S is 
not a group). The first basis is of course {s}sg5, and multiplication in C5' 
with respect to this basis is just the linear extension of the multiplication in 
S. To motivate the second basis, recall that every (finite) inverse semigroup is 
isomorphic to a sub-semigroup of a rook monoid and can therefore be viewed as 
a collection of partial functions. There is another model for composing partial 
functions: only allow the composition if the range of the first function "lines 
up" with the domain of the second. For example, if 

1234\ /1234 
2 - 1 - j' 4 3 - - 

then the idea is that the composition tt o a is 

12 3 4 

3 - 4 - 

and the composition a on is disallowed. The groupoid basis for CS* encodes this. 

Definition (groupoid basis). Let S be an inverse semigroup. Define, for each 
s £ S , the element [sj G CS by 



teS:t<s 

Theorem 3.1. The collection {\_s\}ses is a basis for CS . Multiplication in CS 
relative to this basis is given by the linear extension of 

[sJ [tj = I L'^^J dom(s) = ran(i) 
1 otherwise. 

Furthermore, the change of basis to the {s}s^s basis of CS is given by Mobius 
inversion: 

s= E w- (7) 

teS:t<s 

Proof. This is |31j . Lemma 4.1 and Theorem 4.2, using our convention that 
maps act on the left of sets. □ 

The viewpoint then is that we have two natural bases for CS, the basis 
{s}sgs, and the basis {[sJ}sgs- Note that, if 5 is a group, then s = [s\ £ CS 
for all s e S. 

3.3 Fourier bases for CS 

Let S be a semigroup such that CS is semisimple. This means that CS* decom- 
poses into a direct sum of simple submodules (i.e. left ideals): 



CS = 0Ci,;. 
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Let 3^ be a complete set of inequivalent, irreducible matrix representations 
of CiS*. Then, according to Theorem 12. 5[ 

0p:C5^0A/,^(C) (8) 
pey pey 

is an isomorphism of algebras. 

There is a natural basis for the algebra on the right: the set of matrices in 
the algebra with the property that exactly one entry is 1 (the rest being 0). 
The inverse image of this set is a basis for CS* called the dual matrix coefficient 
basis for y, or the Fourier basis for <CS according to y. When we refer to a 
Fourier basis for CS, we mean any basis of <CS that can arise in this manner by 
choosing an appropriate y. Note that there is a unique Fourier basis for (up 
to ordering) if and only if every irreducible representation of S has dimension 
1 (as is the case for S = Z/nZ, for which the isomorphism dH) is the usual 
discrete Fourier transform and the associated Fourier basis is the usual basis of 
exponential functions) . 

Consider the natural basis for the algebra on the right. The preimage of a 
single column of these elements from the p*'' block is a basis B for a submodule 
of CS, and each element of CS acts on B exactly as described by p. We therefore 
have that i? is a basis for an irreducible submodule of CS* (isomorphic to the 
representation module for p). Since the map above was an isomorphism, the 
preimages of distinct columns have intersection {0}, and we therefore have the 
well-known fact: 

Theorem 3.2. Each irreducible submodule of CS occurs in the decomposition 
of CS into irreducibles exactly as many times as its dimension. 

4 The Fourier transform on CS 

Let S be an inverse semigroup, and let y be any complete set of inequivalent, 
irreducible matrix representations of CS. 

Definition. The isomorphism ([8]) is called a Fourier transform on CS. 

Definition. Let f G CS . The Fourier transform of f is the image of f in the 
matrix algebra 

0Mrf,(C) 
pey 

via the isomorphism ^ . Equivalently, the Fourier transform of f is the re- 
expression of f in CS in terms of a Fourier basis for CS . 

Let / G CS be given with respect to one of the natural bases, i.e., either 
/ = X]s6s/('5)* or / = X^seS L'Sj- We shall sometimes say "calculating 
the Fourier transform on CS"' to mean calculating the Fourier transform of an 
arbitrary element / G CS* given with respect to one of the natural bases, where 
the choice of y is understood. 
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Definition. Let p be a representation ofC-S. Define 




For p G y, f{p) is therefore just the p^^ block in the image of f in the isomor- 
phism ([8]). 

Remark: Since the map in ([8]) is an isomorphism, it respects multiplication, 
and hence the Fourier transform turns convolution of elements of CS* ([T|) into 
multiplication of block-diagonal matrices. It turns convolution into pointwise 
multiplication if and only if all irreducible representations of 5* have degree one, 
which, for example, is the case for S = Z/nZ. 

5 Decomposition of Ci?„ into a matrix algebra 
over group algebras 

Let S be an inverse semigroup. The theorem in this section is a special case of 
Theorem 4.6 in Steinberg [31 j . which provides an explicit isomorphism between 
CS* and a direct sum of matrix algebras over group algebras. The groups ap- 
pearing in the group algebras are the maximal subgroups of S. The purpose of 
this section is to explain this isomorphism in the case when S = 

This turns out to be a very important isomorphism for the FFT theory for 
Rn- It allows us to easily describe a complete set of inequivalent, irreducible 
representations for i?„ (Section l6.ip . it allows us to describe explicit Fourier 
bases for Ci?„ (Section l6.2p . it allows us to easily state the Fourier inversion 
theorem for Ci?„ (Section 16. 3p . and it forms the basis for one of our FFT 
algorithms (Section [9|) . 

Let Dk C Ci?„ be the C-span of {[sj : s G i?„,rk(s) — k}. By ([6]), we have 
that Ci?„ — ©"^0 so that the product of two elements in Dk is an element 
of Dfe, and the product of an element in Di with an element in Dj is if i 7^ j. 

For k G {0, . . . , n}, we identify Sk ^ Rn with the subgroup of elements of 
rank k with domain and range equal to {1, . . . , fc}. We therefore take Sq to be 
the set consisting of the zero map, and = 5*1. The following decomposition 
theorem was implicit in the work of Munn and made explicit in [31j , Theorem 
4.6. 

Theorem 5.1. Dk = Af(^)(C5fc), and thus CR^ = ®l=o M(^,^^{CSk)■ 
Proof. Use the /c-subsets of {1, . . . , n} to index the rows and columns of 

M(„)(C5fc) 

in such a way that {1, . . . ,k} is the first fc-set. If A and B are two fc-subsets of 
{1, . . . ,n}, let P(A~*B) denote the unique order preserving bijection from A to 
B. Keeping in mind that our maps act on the left of sets, define a map 
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by defining it on a basis clement [sj : 

P(ran(s) — . • -^fc}) ^^'(i 1,- ■ -j^} — ^dom(s)) -^ran(s) ,dom(s) ; 

where i?ran(s),dom(s) IS the (^) X matrix with a 1 in the ran(s), dom(s) posi- 
tion and elsewhere. Observe that P(i.an(s)^{i,...,fc})sp({i,...,fc}^dom(s)) e Sk by 
construction. 

It is easy to show that if) is an isomorphism and that is induced by 

gEA^B ^ [P({l,...,k}-,A)9P{B^{l,...,k})\ ■ 

Note that if g G Sk Rn , then 

0(L9J) =.9^1.1 

□ 

Notational remark: For the rest of the paper, if A and B are fc-subsets of 
{1, . . . , n}, then we will denote the unique order preserving bijection from A to 
B by p(^^B). 



6 Consequences of Theorem 15.11 

6.1 Explicit matrix representations of Ci?„ 

Let G be a finite group. The representations of the matrix algebra M„(CG) were 
studied as early as 1942 by A.H. Clifford in the context of Brandt groupoids [3]. 
In the notation of Section O we have 

n 

Ci?„-0Af(„)(C5,). 

fe=0 

Given an irreducible matrix representation p of Sk (or of CSk), we can "tensor 
up" to an irreducible matrix representation p of M^-n-^(CSk) and then extend p 
to an irreducible matrix representation of Ci?„ by declaring it to be on the 
other summands. Specifically, for g G Sk, 

p{gEAM)^ Ea,b® p{g)- (9) 

Let IRR{Sk) be any complete set of inequivalent, irreducible matrix representa- 
tions of Sk- For pi, p2 G IRRiSk), [3] shows that pi and p2 are equivalent if and 
only if pi and p2 are equivalent, and that all irreducible matrix representations 
(up to equivalence) of Mii^{<CSk) are obtained in this manner. 

Therefore, the distinct irreducible representations of Ci?„ are in one-to-one 
correspondence with \S2^q IRR{Sk)- Many explicit descriptions of IRR{Sk) 
are known. Two well-known computationally advantageous ones are Young's 
seminormal and Young's orthogonal forms. A description of the former may be 
found in [2| , and a description of the latter in [9] . 
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To be completely explicit, suppose / G Ci?„, with / = J2seR L*J' ^^'^ 

let p £ IRR{Sk)- Then p is an irreducible matrix representation of Ci?,, , and 

rk(s) = fc 

" E E J2 •l'^Pi{^.■■■.k}^A)9PiB^{l....M))i^A,B(Spig))■ 

AC{l,...,n} SC{l,...,n} gGSfc 
|A|=fc \B\=k 

We have thus effectively described the irreducible representations of Ci?„ by 
describing their actions on the { [sj } basis. We remark that, if desired, one could 
obtain a description of the irreducible representations of i?„ by looking at these 
representations on the {s} basis via ([7|). We also remark that there are other 
ways to describe a complete set of inequivalent, irreducible representations for 
Rn that do not involve decomposing it into a sum of matrix algebras over group 
algebras, or indeed referencing the {[sJ} basis at all. For example, see Grood's 
description in [T^ , Halverson's description in [T3] , or Section 110.21 

6.2 Explicit Fourier bases for Ci?„ 

Let IRR{Sk) be a complete set of inequivalent, irreducible matrix representa- 
tions for Sk, for fc G {0, . . . , n}. For each p £ IRR{Sk), let p denote its extension 
(via the method discussed in Section [6T|l to Ci?„. We take y = [Jp- 

It is now easy to explicitly describe the Fourier basis for Ci?„ according to 
y in terms of the natural { [s\ } basis. That is, if _B C Ci?„ is the set of inverse 
images of the natural basis of ®p^y (C) in the isomorphism 

0p:Ci?„->0A/rf,(C), (10) 

then for each b £ B, 

E ^(^) w- 

sGfl„ 

We will now describe the coefScients b{s). 

We begin by assuming we have an explicit description of a Fourier basis for 
CSk for each k G {0, . . . ,n}. That is, if C is the set of inverse images of the 
natural basis of the algebra on the right in the isomorphism 

p:CSk^ A/,^(C), (11) 

peIB.R{Sk) pelRRiSk) 

then, for each c G C, 

c = E c(a;)a;, 
xeSk 
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and we assume that we know the coefRcients c(x). They may be found, for 
example, by using the standard Fourier inversion theorem for groups (Theorem 

[QD . 

Now, fix p e IRR{Sk). Fix c,j e CSk, 

Cij — y ^ Cij{x)x, 
xGSk 

to be the inverse image in the isomorphism of the element of 

pelRRiSk) 

that is 1 in the i,j position in the p block and elsewhere, p is a block matrix 
whose rows and columns are indexed by the fc-subsets of {1, . . . , n}, and whose 
entries are themselves dp x dp matrices. Call the m*'' fc-set M and the n*^ k-set 
N. 

Theorem 6.1. Let X be a dp x dp matrix with a 1 in the i,j position and 
elsewhere. The inverse image in the isomorphism (jlOp of the element of 
®p^y (C) that is Em.n ® X in the p block and elsewhere is 



Proof. If 7 G iy, 7 ^ p, then 

7 b({i.-,fc}-M) 



J XI "^^J^^) W [P(N-,{l,...,k])\ 



because 



J X! '^^j(^) W \v(N~,{\,....k})\ = 
7( E ^^^■(^) W J -0. 



On the other hand, 

p{\V{{\,...M}^M)\ b(JV^{l,-,fc})J j 

= {E^^x ® Idp) {Ei^i X) {Ei^„ (E) IdJ = E,n.n ® X. 

6.3 The Fourier inversion formula for Ci?„ 
We begin by recalling the Fourier inversion theorem for groups. 



□ 
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Theorem 6.2. Let G be a finite group, and f = X^sgg fi^)^ ^ 

IRR{G) be a complete set of inequivalent, irreducible matrix representations for 

G. Then ^ 

' ' peIR.R(G) 

Proof. See [5S], Section 6.2. □ 

We now state a preliminary lemma, which is a direct consequence of Theorem 
[Oand Section O 

Lemma 6.3. Let f — X^xei? /(^) L^^J ^ CRn, and let y be a complete set of 
inequivalent, irreducible matrix representations for i?„, induced by 

l+J IRRiSk) 

ke{0,...,n} 

in the manner described in Section l6.ll Let rk(a;) = k. Write 

X = P({l,...,k}~,ia.n(x))yP(doi-a{x)~,{l,...,k}) 

for a unique y € Sk- 

The ran(a;), dom(a;) entry of f{p) is a dp x dp matrix. Lf we denote it by 
/(p)ran(x),dom(2:), then wc havc 

f{x) = X! dptVACe (/(p)ran(:c),dom(2;)P(y"^)) • 

Proof. For p e IRR{Sk) we have 

with p(s) = if rk(s) ^ fc. By definition of p, the Fom'ier transform /(p) may 
be computed block by block (and thus may be inverted block by block). The 
ran(x), dom(a;) entry oi f{p) is determined by the f{s) for which ran(s) = ran(x) 
and dom(s) — dom(x), and such /(s) do not affect any other entries of f{p). 
Explicitly, the ran(a;), dom(a;) entry of f{p) is given by 

f{P{{l,...,k}^Tan(xy)SP{dom{x)^{l,...,k}))pis). 

seSk 

Let us define a function fran{x),dom{x) on Sk by 

/ran(a:),dom(a;)(s) = f {P({1,. . .,k}^ian(x))SP(dom{x)~>{l,. . .,k})) ■ 

Then 

E fiPi{l,---,k}^ran(x))SP(dom{x)-.{l,...,k}))p{s) = ^ fran{x) Aom{x){s) p{s) 
sGSfc seS'k 

/ran(a:) ,dom(2:) 
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The Fourier inversion theorem for groups then appHes, and yields: 

fi^) = /(P({l,...,fc}^ran(2:))J/P(dom(a:) — {l,...,fc})) 

^ rfptrace (^/ran(2;),dom(a;)(p)p(y~^)) , 

and since 

yran(a;) ,dom(a:) (P) = /(p)ran(a;),dom(a:), 

we are done. 



1 

I'S'fcl 



□ 



Now, let S be any set of inequivalent, irreducible matrix representations for 
Rn- We define some notation before stating the Fourier inversion theorem for 

Let 3^ be a complete set of inequivalent, irreducible matrix representations 
for Rn induced by l+J^j-g^Q „} IRR{Sk) in the manner described in Section [01 
If p G IRR{Sk), we have the corresponding p € y, which is equivalent to some 
representation in which we denote by p. 

Theorem 6.4 (Fourier inversion theorem for Ci?„). Let 



Let rk(x) ~ k, and let us denote the semigroup inverse of x by . Then 

otrace(/(p)p([2:"^J)) . 

Proof. Since p is equivalent to p, write 



' pGlRR(Sk 



p^A-^pA 

for some invertible matrix A. We therefore have 

f{p) = A-^f{p)A. 
As in Lemma 16.31 let y be the unique element of Sk such that 

X = P({l,...,fc}^ran(2:))2/P(dom(x)^{l,...,fc})- 

Now, we have 



trace 



(/(p)ran(x),dom(3;)P(2/ ^! 



=trace(^ /(p) 
=trace (/(p)p([x"^J 
=trace 
=trace 
=trace 



dom(x) ,ran(2:) 



'p(y-')]) 



A-'f{p)A\ [A-'p{[x-'\)A] 
(A-^f{p)-p{[x-'\)A 

(/(p)p(L^-'J 



the last equality arising from the similarity-invariance of trace. The theorem 
now follows from Lemma [631 D 
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7 Associating functions on S to elements of CS 



We are now ready to address the issues that arise in choosing how to associate 
functions on S to elements of the semigroup algebra CS. We also define the 
Fourier transform and the convolution of functions on S. 

Let S be an inverse semigroup. Let / : 5 C. Having defined two "natural" 
bases for CS, we have two natural choices for how to associate / to an element 
of CS, either 

f^J2f{s)s or f^J2f{s)[s\. 
ses ses 
In the first case, which we shall call the semigroup association model, the 
elements s in CiS* are associated to the characteristic functions of the elements 
s G S*. In the second case, which we shall call the groupoid association model, the 
elements [s\ in CS are associated to the characteristic functions of the elements 
s€S. 

Definition (Fourier transform of a function on S). Let f : S ^ C. Given 
an association model, the Fourier transform of f is defined to be the Fourier 
transform of the associated element of CS . 

Thus, if the groupoid association model is used, then the Fourier transform 
of / is 

peyses 

and if the semigroup association model is used, then the Fourier transform of / 
is 

peyses 

where 3^ is a complete set of inequivalent, irreducible matrix representations of 
CS*. Convolution of functions on S is also defined by the association model. 

Definition (convolution of functions on S). Let f,g : S ^ C. Choosing an 
association model defines how the images of f and g in CS multiply. Denote 
the images of f,g in CS by f,g respectively. Then we define, for s G 5; 

Under the semigroup association model, f *g{s) — the s*'' coefficient of f * g 
expressed with respect to the {s} basis. Under the groupoid association model, 
f * g{s) = the [s\*^ coefficient of f*g expressed with respect to the {[s\} basis. 

For the case S ~ Rn, what follows is an overview of the considerations one 
should take into account when choosing which function association model to 
use. 

Consideration 1 (Convolution). Convolution of functions on i?„ is defined by 
the association model used. If the semigroup association model is used, then 
convolution of C- valued functions /, g on _R„ is defined by U) . Specifically, 

{f*9){s)^ f{r)g{t). 

r,t^Rn:rt—s 
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If the groupoid association model is used, then convolution of C-valued func- 
tions /, g on Rn is defined by their multiplication in Ci?„ with respect to the 
{[sj} basis. Specifically, 

ran(s)— raii(r) 

Consideration 2 (Inner products). There is a natural inner product for C- 
valued functions f,g on i?„: 

It would be nice to use this inner product. Choosing an inner product on 
functions and an association model induces an inner product on Ci?„ (or equiv- 
alently, choosing an inner product on CRn and an association model induces 
an inner product on functions on i?„). Specifically, if we choose the above in- 
ner product, then the semigroup association model induces the inner product 
< •, • >i on CRn defined by 



< s,t >i = 



1 if .s = i 
otherwise, 



and the groupoid association model induces the inner product < ■, ■ >2 on Ci?„ 
defined by 

'l if ls\ = lt\ 



< W , W >2 = 



otherwise. 



However, choosing the inner product < ■, ■ >i leads to undesirable consequences 
with regards to the Fourier basis for Ci?„, as described in Consideration [3l For 
this reason, if the semigroup association model is used, then it is recommended 
to take the inner product < •,• >2 on CRn and use that to induce an inner 
product on C-valued functions on i?„. 

Consideration 3 (Orthogonality of isotypic subspaces of Ci?„). Since CRn is 
semisimple, we have 

CRn = Li, 

where each Li is simple. Let us group them according to their isomorphism 
classes: 

CRn = Vp, 

where 3^ is a complete set of inequivalent, irreducible matrix representations 
for Rn, and Vp is the sum of all simple submodules of Ci?„ isomorphic to the 
representation module for p e 3^. The Vp are called the isotypic subspaces or 
isotypic submodules of CRn- Note that this decomposition does not depend on 



21 



the choice of 3^. Let v £ Vp,v' £ Vp'. We would Uke to have an mner product 
(•, •) on CRn such that, if Vp ^ Vp' , then 

i.e., an inner product under which the isotypic subspaces of Ci?„ are orthogonal. 
The claim is that < ■, ■ >2 accomplishes this, and that < ■, ■ >i in general does 
not. 

Theorem 7.1. Let Vp ^ Vp> . Then, using the notation above, < v,v' >2— 0. 

Proof. By linearity, it suffices to show this in the case that v and v' are Fourier 
basis elements of Ci?„. We take y to be the set induced by l+)feg|o «} IRR{Sk) 
in the manner described in Section l6.ll and we assume that v and v' are part 
of a Fourier basis for Ci?„ according to y. 
We know that 

CSk = Wp, 

pelRRiSk) 

where Wp contains all the irreducible submodules of CSk isomorphic to the 
representation module for p. If w £ Wp,w' € Wp',Wp ^ Wp', then imder the 
inner product [•, •] on CSk defined by 

it follows from the discussion in Chapter 2 of 28J that we have [w, w'] = 0. 

Now, suppose that p was induced by p G IRR{Sk) and that p' was induced 
by p' £ IRR{Sj). By Theorem 16. 1[ when written in terms of the natural {[sj} 
basis, V contains nonzero coefficients only for the elements \_s\ where s is of rank 
k, and v' contains nonzero coefficients only for the elements [sJ where s is of 
rank j. Thus, if A: ^ j, we have < v,v' >2— 0. Suppose then that k — j. By 
Theorem 16. H we have 

V = b({l,...,fc}^B)J W [PiA^{l,...,k})\ ' 

s&Sk 

v' = [P{{U...,k}^B')\ ^'(^) W [PiA'^{U...,k})\ ' 



for A, A' , B, B' some fc-subsets of {1, ... , and 

Y <s)s £ Wp, Y ^'(s)^ e Wp' 
seSk ssSfc 

some Fourier basis elements for CSk- 

li A ^ A' or B ^ B' , it is apparent that < v,v' >2— 0, so suppose further 
that A = A' and B = B' . 
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Now, since Vp ^Vp', we have p^^ p', and we therefore note that 

[J2v{s)s,J2v'{8)s]=0. 

seSk seSk 

For convenience, we temporarily denote P({i,...,fe}^B) by pb and P(A->{i,...,fe}) 
by Pa- We have 

< v,v' >2= ^ ^ ■u(s)t;'(t) < [pbspa} , bsip^J >2, 

seSfc teSfc 

and [pssp^J = [PstpAj if and only if s = t, so 

< V, v' >2= ^ v{s)V{s} = [ J2 v{s), J2 ^'(«)] = 0- 

seSk seSk seSk 

□ 

Remark 7.2. The isotypic subspaccs Vp need not be orthogonal under tlic inner 
product < •, • >i. For example, consider CRi, which has two nonisomorphic 
irreducible representations, each of degree 1, and hence a unique Fourier basis. 
It decomposes as 

CRi = (C - span([ldj)) (C - span([Oj)) . 

We have 

< LidJ , LOJ >!=< Id - (0), (0) >i= -1. 

8 FFT algorithm considerations 
8.1 Computational complexity 

If G is a finite group, then {g}g^G indexes the natural basis of CG. If 5" is an 

inverse semigroup, then CS has two natural bases, the {.s},5gs basis and the 
{ [^1 }seS basis. We therefore define the following two notions of computational 
complexity for the Fourier transform on CS. 

Definition (Computational complexity). Let y be a complete set of inequiva- 
lent, irreducible matrix representations forCS. For an arbitrary element f G CS 
expressed with respect to the {s} basis 

the minimal number of operations to compute the Fourier transform of f, i.e., 
to compute 

f{p) = Y,fis)p{s) (12) 

ses 
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for all p £y, is denoted by Ty{S). 

For an arbitrary element f G <CS expressed with respect to the { [sj } basis 

seS 

the minimal number of operations to compute the Fourier transform of f , i.e., 
to compute 

for all p € y, is denoted by 7^"^ (S). 

Now, let y vary over all complete sets of inequivalent, irreducible represen- 
tations for <CS . We define 

C%S) = inmy{T^{S)), 

CW (S) = miny{Tj'^ (S*)), and 
0(5) =min(C"(5'),CW(S')). 

An operation is defined to be a single complex multiplication followed by a 
complex addition. For the purposes of our analysis, we assume that all repre- 
sentations in y are precomputed and stored in memory. 

If G is a group, the {s} and the { [sJ } bases of CG are identical, so we may 
drop the superscripts on the complexity notation. For example, 

Now, let y be any complete set of inequivalent, irreducible matrix represen- 
tations for CRn. A naive implementation of the Fourier transform on i?„, i.e., 
computing and directly, gives 

C{Rn) < El^»I^P' 

which, by gives 
Theorem. C(i?„) < |i?„p. 

As n grows, this cost quickly becomes prohibitive. 

For the same reason, we also have C(G) < |Gp for any group G, but many 
families of groups enjoy results along the lines of C(G) = 0(|G| log'^ |G|). Indeed, 
such upper bounds remain the goal in group FFT theory. There are also groups 
G for which there currently exist greatly improved (but not 0(|G| log'^ |G|)) 
algorithms, such as matrix groups over finite fields or, more generally, finite 
groups of Lie type |27| . It is conjectured that there are universal constants 
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Ci,C2 such that for any group G, C{G) < ci|G[log°^ \G\ "51]. It is currently 
known that, for example, 

C(Z/nZ) = 0(n log n), 

C(G) < 8.5|G| log |G| for any supersolvable group G, and 

C{Sn) = 0(|5'nllog^ I'S'nl) (wherc 5„ is the symmetric group on n letters). 

See [B], [T], and [TH], respectively, for these results, or see Section 1531 for a 
description of an 0(|5„| log^ |S'„|) FFT on S„. 
We now state the main results of this paper. 

Theorem 8.1. CW(i?„) < 0(|i?„| log^ \R„\). 

Theorem 8.2. For any e > 0,C''(i?„) = 0{\Rn\'^+'), and hence 

C{Rn)^0{\Rn\^+'). 

We prove these results in Sections 19.11 and 19. 2i respectively. 

8.2 FFTs on groups 

In this section, we describe a general method used to construct FFTs on groups. 
This section is basically an elaboration of the main idea in [37] , and is included 
in the interest of keeping this document self-contained. This method is used in 
the FFT on S'„ given in Section [5751 (which itself is used in the FFT for i?„ given 
in Section [9]). We also will generalize this method directly to create an FFT for 
Rn in Section [TOl 

Let G be a group and let _ff be a subgroup of G. H partitions G into cosets 
yH. Let G C G be a full set of coset representatives for distinct cosets (i.e., 
G = Uy6C ^'^'i 2/1 7^ 2/2 e G yiff n 2/2 -ff = 0)- Let p be a representation 
of G. Then, since p is a homomorphism, we have the factorization 

fip) - E - E p^y) E fy^h)p^h), (14) 

gSG yeC h&H 

where fy{h) = }{yh). 

Let ya and yu be complete sets of irreducible, inequivalent matrix represen- 
tations for G and i?, respectively, adapted to the chain G > H. We would like 
to compute f{p) for all p e yc- The idea is that if we already knew all the fyi^j) 
for all 7 in yn and y G C, then we could construct the inner sums in (fT4)) based 
on how p splits when restricted to H. (To do this in practice for a particular 
group G, we would need a theorem that says how p splits when restricted to H. 
Such theorems are known as branching theorems.) Since yn, are adapted, 
this construction can be done for free, since (say p\h — pi (B ■ ■ ■ (B Pk) 

E fv^^)P^^) = E ^y^^^ [^i(^) ® ■ ■ ■ ® ^fe(^)] = ® • • • ® fyiPk)- 

heH heH 
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Then, to finish computing f{p), we need to multiply matrices together and 

add the results. Denote J2heH fyWpi^) ^a(p)- Then by choosing our 
representations yc, for G and H to be adapted to the chain G > H, we 
have 

Theorem 8.3. 

TyJG)<^-^^Ty„{H) + MG{C), 
where Mq{C) is the number of operations needed to compute the sums 

y&C 

for all p G yc, given the matrices Ay{p). 

Mg{C) can be made small by choosing the coset representatives y £ C in 
such a way that p{y) is a sparse matrix (or, more generally, can be factored into 
a small number of sparse matrices). 

Of course, this method can be used recursively with a subgroup chain and a 
corresponding collection of chain-adapted representations, as in [^T], [57]. 

8.3 An FFT on the Symmetric Group 

In this section, we describe an algorithm to compute the FFT on Sn- The FFT 
algorithm for Rn presented in Section [9] exphcitly uses FFTs on Sn, and ideas 
from this algorithm are used in the FFT on i?„ presented in section [TOl 

Consider Sn and the subgroup chain Sn > Sn-i > ■ ■ ■ > Si — {!}, where 
Sk is identified with the subgroup of Sn that fixes the points k + 1, . . . ,n. 
Choose a complete, inequivalent, irreducible set of matrix representations for 
CSn adapted to this chain and call it 3^„. Two well-known choices are Young's 
orthogonal and Young's seminormal forms (see, for example, [S]). Let tj be the 
transposition {j — Let e denote the identity element of Sn- We use the 

following set of coset representatives: 

{Ti : 1 < i < n,Ti = e_^^^ti+iti+2 ■ ■ ■ tn}- 

i 

Thus, for example, Ti = ^2^3 ■ ■ - tn, ^n-i — tm and T!„ = e. We therefore obtain 

n 

J2f{s)p{s)^J2p('^^^ E fTMpis), (15) 

which then implies 

TyJSn) < -^ry_,{Sn-l) + A^S„({T,}). (16) 

Once we have calculated all the /t; on Sn-i, then for any p G yn, we can 
reconstruct (for free, since we are using chain-adapted matrix representations) 
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the inner sums in based on how p sphts when restricted to Sn-i- This 
sphtting is described by the branching theorem for S'„ (see, for example, [9], p. 
304). 

We now turn to analyzing the Ms^({Tj}) term in (fTB|). Notice that, for 
j > 2, tj e Sj and tj commutes with Sj-2 (and t2 £ S2 and ^2 commutes 
with 5*1). It is easy to derive from the combinatorics of Young tableaux (see, 
for example, [15j ) that the maximum multiplicity occurring in the restriction of 
any irreducible representation of Sj to Sj-2 is 2. By Schur's Lemma, then, for 
any j (2 < j < n) and any p G 3^„, p{tj) contains at most 2 nonzero entries per 
row and column. 

Fix p G 3^„. Since p{Ti) = p{ti+i) p{ti+2) ■ ■ ■ p{tn) , computing p{Ti)ATXp) 
for an arbitrary matrix At^ (p) may be accomplished by multiplying Ay. {p) on 
the left by p(tn), multiplying the result on the left by p(tn-i), multiplying 
the result of that on the left by p{tn-2), etc. It therefore takes a maximum 
of 2(n — operations to perform the multiplication p{Ti)ATi{p) (keeping in 
mind that multiplying by the identity matrix may be done for free), and once 
p{Ti)ATi{p) has been computed for all Ti, it takes a maximum of (n — l)(i^ 
operations to add the results to give jy^=i P(^i)^Ti (p)- Letting p vary over 3^„, 
then, implies 

n 

Ms Aim < J2Y.'^{n~i)dl+Y,{n~l)dl 

= {n~\){n)\S^\ + {n^\)\S,,\ 
= (n2-l)|5„|. 

Putting this together with (fT6|) . we obtain: 

Ty^{Sn) < nTy^_,{Sn-i) + {n^ - l)\Sn\. 
Induction on n then yields 
Theorem 8.4. 7^„(S'„) < |n(n + lfn\ 

This is of order n!(logn!)'^ = \Sn\ log'^ |S'„|. The algorithm described here is 
the heart of Clausen's FFT on 5„ [J. 

By a more careful analysis of the matrix multiplications involved, Maslen 
has obtained an algorithm for the FFT on Sn of complexity 0(|5'„|log^ l^nl) 
([18], Theorem 1.1): 

Theorem 8.5. Let denote a complete set of inequivalent, irreducible repre- 
sentations for Sn in Young 's orthogonal or seminormal form. Then 

TyASn) < ^n{n-l)\Sn\. 

9 An algorithm for the FFT on 

In this section, we present the faster of our two algorithms for the FFT on i?„. 
In the process, we prove Theorems 18.11 and 18.21 
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9.1 An FFT for changing from the {[sj} basis to a Fourier 
basis 

Let / e Ci?„ be an arbitrary element, given with respect to the {[sJ} basis: 

/ = E /(^) w ■ 

Let yk denote Young's seminormal (or orthogonal) matrix representations 
for Sk. For our complete set of inequivalent, irreducible representations for i?„, 
we take the set y induced by l+Jfcg{o n} -^fc manner described in Section 

16.11 U p G y was induced hy p £ yk, then we have that f{p) is an (^) x (^) 
matrix whose rows and columns are indexed by the fc-subsets of {1, . . . , n}, with 
entries themselves dp x dp matrices. By Theorem 15.11 and we know that for 

s e Rn, 



p(W) 



if rk(s) ^ k 

^Eran{s),dom{s) ® P(P(ran(s)-.{l, . . . ,fc}) SP({1, . . .,fc}-,dom(s)) ) Otherwise. 

Let A and B be fc-subsets of {1, . . . , n). Then the A, B entry of f{p) is 

f{p)A,B = E fiP({^,---M^AsPB^{i,...,k}))P'y^)- 

seSk 

If we define a function on Sk by 

fA.nis) = fiP{{l,...,k}~>ASPB^{l,...M}))' 

then p7p is just fA,Bip), a Fourier transform on 5*^. An obvious algorithm 
presents itself: for each k, run (^)^ FFTs on 5*^. By Theore: 

Ty,{Sk)<lk{k~l)\Sk\. 

Theorem (Theorem [H]). CW(i?„) < 0(|i?„|log^ |i?„|). 
Proof. By the above algorithm. 



n / \ 2 

k=0 



k, 



k=0 

< ^n{n-l)\R^\. 



Since (i?„) < T-^'i > n\, and n = 0(log(n!)), we are done. □ 
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9.2 An FFT for changing from the {s} basis to a Fourier 
basis 

Let / e Ci?„ be an arbitrary element, given with respect to the {s} basis: 

/ = E f('>- 

sGfl„ 

Thanks to the algorithm in Section 19.11 a simple algorithm for calculating 
the Fourier transform of / is evident: use the set y of incquivalent, irreducible 
matrix representations for i?„ used in Section 19.11 re-express / in terms of the 
{ [sj } basis, and then run the algorithm in Section 19.11 Under this approach, 
the only algorithmic complexity left to consider is the complexity of changing 
from the {s} basis to the {[sJ} basis. Suppose 

/= E 9i-s)ls\. 



Since 



basic linear algebra gives 



= E 

teR„:t<s 



gis) = E /(^) 



Given s G Rn with rk(s) = k, there are \Rn-k\ elements a: G i?„ satisfying 
X > s. A naive implementation of the {s} — > {[sJ} change of basis therefore 
takes at most 



E 

fe=0 



2 

n 



kl\Rn-k\ (18) 



operations. 

If we consider the matrix that encodes this change of basis (which consists 
only of O's and I's), then (fT8|) is the number of I's in this matrix, counted row 
by row, since a row corresponding to an element s S Rn of rank k has |i?„_fc| 
I's in it. On the other hand, a column corresponding to an element s G i?„ of 
rank k has 2*^ I's in it, so counting column by column, we obtain: 

n / s 2 Ji/\2 



A;!|i?„_fe| - E J ^'2 



k=0 ^'^^ fc=0 

and thus the {s} change of basis takes at most 

EU) fc!2'=<2"EU) 

fc=0 ^ ^ k=0 ^ ^ 

operations. 

Putting this together with Theorem 18.11 we obtain 
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Theorem (Theorem [OI . Lete>0. Then 

Proof. The algorithm in this section consists of two steps: change from the {s} 
basis to the {[sj} basis, and then rmi the algorithm in Section [?7T] Thus 

<2"|i?„| + r3;W(i?„) 

<2"|i?„| + ^n(n-l)|i?„|. 

For large enough n, we have 2" > — 1) and > n!*^ > 2", so 

T^{Rn)<0{\R,,\'+'). 

□ 

Remark: We have developed improved algorithms for the {s} {[s\} change 
of basis, but none yet which are 0(|i?„| log^ \Rn\)- Also, in Section [10] we exhibit 
a different algorithm for changing from the {s} basis to a Fourier basis, which 
avoids the {[sJ} basis entirely. However, it too is a 0(|i?„p+^) algorithm. 

10 Another algorithm for the FFT on Rn 

In this section, we present an algorithm for the FFT on i?„ that is quite dif- 
ferent from the algorithm presented in Section [31 The algorithm in Section [5] 
involved examining the poset structure of i?„, decomposing Ci?„ into matrix 
algebras over group algebras, using matrix representations for i?„ induced by 
computationally advantageous matrix representations of Sk , and running FFTs 
on Sk- 
in contrast, the algorithm in this section directly generalizes the ideas in 
Sections 18.21 and 18.31 It involves finding a semigroup factorization of i?„ into 
"cosets" and using a chain-adapted set of matrix representations to run recursive 
calculations on the sub-semigroup chain i?„ > Rn-i > ... > Ri- The algorithm 
in this section computes the Fourier transform of an element in Ci?„ expressed 
with respect to the {s} basis directly, without referencing the {[sJ} basis. It 
is not an 0(|i?„| log'^ |i?„|) algorithm, so we do not consider the {[sJ} {s} 
change of basis necessary to run it on the { [sJ } basis. Even for an element 
expressed with respect to the {s} basis, the algorithm in Section [9] is slightly 
better. Nevertheless, we present this algorithm because the ideas involved may 
be helpful for designing FFTs on other inverse semigroups. 
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10.1 Another FFT for changing from the {s} basis to a 
Fourier basis 

Let / e Ci?„ be an arbitrary element, given with respect to the {s} basis: 



As with the symmetric group (Section 18. 3p . let tj be the transposition 
(j — 1, j). We use the following sets of "coset representatives" (we call them 
"coset representatives" because they play the same role in this FFT as coset 
representatives do in group FFTs): 

{Ti : 1 < i < n,Ti — U+iti+2 ■ ■ ■ tn} (where T„ — Id) and 

{TM < i < n - 1, r = tntn-l ■ ■ ■ k+i}. 
We use the semigroup chain i?„ > i?„_i > . . . > where 

-Rfc = {cr G i?„ : a{j) = j if j > k}. 

Halverson |13| has found a complete set of inequivalent, irreducible matrix rep- 
resentations for Rn adapted to this chain. The description may be found in 
Section flO. 21 Call this set y„. For each p G 3^„, we must compute 



A subsemigroup does not necessarily partition its parent semigroup into 
equally sized cosets, so we cannot directly factor through a subsemigroup as in 
(jl4p . Instead, we use an approach for R„ that is based on the recursive formula 
([2]). With this, we have the following factorization theorem. 

Theorem 10.1 (Factorization theorem for For any representation p of 

Rn, if n > 3, we have the following factorization. 



f= E /(^)^- 



n 



f{p) = Y^p{T,) E /T.(s)p(s)+p([n]) E /m('^)p(^) (19) 




ri-l 



+ E E f\')p(') p(T'), 



where [n] is the link (1)(2) . . . (n — ^)[n], /yi(s) — f{As), and 




if n — 1 G dom(s) 

/(sT*) otherwise. 



Proof See Section [inSl 



□ 
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As in the group case, we will use the above breakdown to compute /(p) 
recursively. If we knew /[„](7) for all 7 £ 3^„-^i, friil) for all 7 G 3^„-i and 

1 < i < n, and /^'(7) for all 7 G 3^„-i and 1 < i < n — 1, using (fT9)) we 
could reassemble them for free, since we are using chain-adapted representations, 
based on how p splits when restricted to Rn-i (see Theorem 110. 6|) to calculate 
f{p) for any p ^ Therefore, we have 

Lemma. For n > 3, 

where Mr^ is the total number of operations required to compute the sum (|19p 
for all p G given knowledge of the /[„], all the fr^, and all the on Rn-i- 

We now analyze Mr^ to obtain: 

Theorem 10.2. For n > 3, 

T^^iRn) < 2nry^_^(i?„„i) + 2n^\R„\. (20) 
Proof. To analyze Mr^, let: 

• Ml = The maximum number of operations necessary to calculate the 
matrix product p{Ti)ATi{p) for arbitrary matrices At^Ip), for all p £ yn 
and for aU (1 < i < n). 

• M2 ~ The maximum number of operations necessary to calculate the 
matrix product A'^ {p)p{T^) for arbitrary matrices A^ (p), for all p G 3^„ 
and for aU (1 < i < n - 1). 

• M3 — The maximum number of operations necessary to calculate the 
matrix product p{[n]) A[n]ip) for arbitrary matrices A[„](p), for all p G 3^„. 

• — The maximum number of operations necessary to add together 2n 
dp X dp arbitrary matrices, for all p G yn- 

Then Mr^ < ^ti M^. 

Analysis of Mi: For each p G yn and each T^, we must perform the multi- 
plication p{Ti)ATi{p) — p{ti^i)p{tij^2) ■ ■ ■ p{tn)ATi{p) for an arbitrary matrix 
ATi{p). As was the case with Sn, for j > 2, tj G Rj, tj commutes with Rj-2, 
and M{Rj, Rj-2) = 2. By Schur's Lemma, p{tj) {j > 2) contains at most 2 
non-zero entries per row and column. For j ~ 2, t2 G i?2, and the maximum 
dimension of an irreducible representation of R2 is 2. Therefore, p(t2) contains 
at most 2 non-zero entries per row and column. Alternatively, it is obvious from 
the description of ^„ (see Section [10.2p that p{tj) {j > 2) contains at most 2 
nonzero entries per row and column. Therefore, multiplying an arbitrary ma- 
trix by p{tj) on the left requires at most 2d^ operations, and so performing the 
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multiplication p{Ti)ATi{p) requires at most 2{n — i)d'^ operations. Therefore, 
we have 

n 

Ml < E E " H = - 

where the final equality comes from ([?]). 

Analysis of M2: The only difference between Mi and M2 is that M2 involves 
multiplying arbitrary matrices by p{T^) on the right rather than by p{Ti) on 
the left, so M2 is the same as Mi in the complexity analysis. Thus 

M2 < [n)[n - l)\R^\. 

Analysis of M3: Since [n] S i?„, [n] commutes with Rn~i, and M{Rn, Rn-i) = 
1, we have that contains at most 1 non-zero entry per row. Thus 

Analysis of M4: For a particular p, the matrix additions can be accomplished 
with (2n — operations. Thus 

Mi < E (2?^ - = i^n - l)|i?„|. 

Putting this all together, we obtain 

Mr„ < {2{n){n - 1) + 1 + 2n - l)|i?„| = 2n^\Rn\. 

□ 

We now prove that this algorithm gives the complexity result 
Theorem 10.3. For n > 5, T^„(i?„) < 2"n|i?„|. 

Proof. Base case: |i?2| = 7, so a naive implementation of the FFT on R2 gives 
TyAR^) < 49. Applying repeatedly, we have 

TyARs) < 2{i)Ty,{R2) + 2{if\R^\ < 6(49) + 18(34) = 906, so 
Ty^Ri) < 2{4)Ty^{R3) + 2{Af\Ri\ < 8(906) + 32(209) = 13936, so 
Ty.iRs) < 2{b)Ty^{Ri) + 2(5)2|i?5| < 10(13936) + 50(1546) = 216660, and 
216660 < 25(5)|i?5| - 247360. 

This proves the base case. 
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Similarly, for rt = 6, we have 

Ty.iRe) < 2{6)Ty,{R5) + 2{6f\R6\ < 12(216660) + 72(13327) 
= 3559464, and 3559464 < 2^{6)\R6\ = 5117568. 

Now, let n > 7. Observe that, for a = 1 to n, the sets {a E Rn : (T{a) = n} 
are disjoint, and each are of size thus n|i?„_i| < Therefore, we 

have 

TyJRn) < 2nTy^_,iRn^i) + 2n^\R„\ 

< 2n(T-^{n - l)|i?„-i|) + 2n^\R^\ 

< 2"(n- l)|i?„| +2n2|i?„| 
= 2"n|i?„| + (2n2_2")|i?„| 

< 2"n|i?„| (since 2n^ < T for n > 7). 

□ 

10.2 Chain-adapted matrix representations for _R„ 

In this section, we give a description of a complete set of irreducible, incquiv- 
alent, chain- adapted matrix representations for i?„ relative to the chain i?„ > 
Rn-i > ■ ■ ■ > Ri- The results herein are a special case of the results in [T5] . 

Definition (partition). A partition X of a nonnegative integer k {written A h fc) 
is a weakly decreasing sequence of nonnegative integers whose sum is k. We 
consider two partitions to he equal if and only if they only differ by the number 
of 's they contain, and we identify a partition A with its Young diagram. 

For example, A — (5, 5, 3, 1) is a partition of 14, and 
A= (5,5,3,1)- (5,5,3,1,0) = 



It is well-known that a complete set of inequivalent, irreducible representa- 
tions for Rn is indexed by the set of all partitions of the integers {0, 1, . . .«} 
(see, for example, [12] or [29]). Therefore, for integers n > 0, let 

A„ {A h fc : < fc < n}. 

Definition (n-tableau, n-standard tableau). For A e A„, define L to be an 
n-tableau of shape A if it is a filling of the diagram for A with numbers from 
{1, 2, . . . n} such that each number in L appears at most once. L is an n-standard 
tableau if, furthermore, the entries in each column of L increase from top to 
bottom and the entries in each row of L increase from left to right. 
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Fix A. Let denote the set of n-standard tableaux of shape A. The 
symmetric group acts on tableaux by permuting their entries. If L is an n- 
tableau, then [i — 1, i)L is the tableau that is obtained from L by replacing i — 1 
(if i — 1 G L) with i, and by replacing i (if i G L) with i — Note that L G 
need not imply {i — 1, i)L e T^. 

Let {vl : i S T^} be a set of independent vectors. We form 

= C-span{wL : L G T^}. 

As such, the symbols Vl, for L G T^, are a basis for the vector space . 
Halverson defines an action of on in such a way that (extending by 
linearity) is an irreducible Ci?,„-mod\ile and such that, as A ranges over 
A„, the V'^ constitute a complete set of inequivalent, irreducible representation 
modules for We first describe this action, and we then describe an ordering 
of the bases for the so that the resulting matrix representations are chain- 
adapted to Rn > Rn-l > ■ ■ ■ > Rl- 

Definition (content). If b is a box of X in position {i,j), then the content ofb 
is defined to be 

ct{b) =j-i- 

Let L G T^. If i — l,i G L, then let L{i — 1) and L{i) denote the box in L 

containing i — 1 and ?', respectively. 

To define the action of i?„ on , it is sufficient to define the action of a set 
of generators of i?„ on 

Definition (action of Rn on V^). Define the action of the transpositions ti = 
{i — l,i), for 2 <i <n, as follows: 



tiVL 



where 



ctiL(i))-ct(L(^-l)) ^L + (1 + ct{L(^))-ct(L{^-l)) >L' if i - l,i & L 

vtiL if exactly one of 

i — l,i E L 

^Vl if i — l,i ^ L 



VL' = 



VtiL if tiL is n-standard 
otherwise. 



Define the action of the link (1)(2) . . . {n — l)[n] = [n] on F'^ by 



[nlVL 



Vl ifn^L 
if n e L. 



Remark: If A = (0), then is 1-dimensional, and the action of i?„ on is 
the trivial action given hy xv = v for all x G Rn and all v £V^. 
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Theorem 10.4. As A varies over all partitions of all nonnegative integers less 
than or equal to n, the constitute a complete set of irreducible, pairwise 
non-isomorphic representation modules for Rn 

Definition (corner of a partition). A corner is a box c of X for which X contains 
no box to the right or below c. In other words, the corners are the possible 
positions of n in an n-standard tableau of shape X. 

We now record the branching theorem for i?„. 

Theorem 10.5 (Branching theorem). As a Ci?„_i module, 

where X~'^ is the set of all partitions /i G A„_i such that either ^, ~ X {if XV n) 
or fi is obtained by removing a corner from X ^ISf . 

Now, for purposes of chain-adaptation, we order the basis {vl} for using 
the following generalized last-letter ordering. 

We begin by partitioning the vl into subsets based on the corners ci , . . . , c; 
of A. Number the corners from top to bottom. Now, form the sets 

V^{0) ^{vL-.Le and n ^ L}, 

V^{i) ^{vl: LeT^ and n G q of L}, l<i<l, 

and declare elements of V'^{j) to be earlier in the ordering than elements of 
V^{k) whenever j < k. To order the subset V^{k), delete the corner Ck (do 
nothing if fc = 0) and repeat the same ordering process (starting by identifying 
the corners of the resulting partition and partitioning the into subsets based 
on those corners) with n — 1 in place of n, etc. 

As an example, consider A = (2, 1, 1), which has two corners, and R^. Our 
ordered basis for the 15-dimensional is 



< ^ 



< ^ 



1 


CM 


3 




4 




1 


4 


2 




5 





1 


5 


2 




3 




1 


4 


3 




5 





< vr 



< vr 



< vr 



< vr 



Remark: If A h n, then the generalized last-letter ordering scheme given above 
reduces to the usual last-letter ordering scheme used for Young's orthogonal and 
seminormal representations of the symmetric group. 

It is now easy to see, under this ordering of the bases for the , that the 
matrix representations described in this section are chain-adapted to the chain 
Rn > Rn~i > ■ ■ ■ > Ri- We know, by the branching theorem for i?„, how 1/'^ 
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decomposes as a Rn-i module (it decomposes into the C-span of the V^{k) for 
< fc < /), and it is obvious by the action of Rn~i on F'*' that 

i?„_iC-span(y^(fc)) C C-span(y^(fc)) 

for all k. The same argument is used to induct down the chain R„ > Rn-i > 
... > Ri, now, and is trivial because we ordered our basis for V'^ inductively 
according to the same rule. We now restate the branching theorem for i?„ under 
our ordering of the bases for the T^'^. 

Theorem 10.6 (Branching theorem). Let be the matrix representation as- 
sociated to with respect to the basis {v^}, with the basis ordered according 
to the generalized last-letter ordering. Then 

where A~'^ is the set of all partitions fi G A„_i such that either p = X or fj, 
is obtained by removing a corner from A. The first fi £ A^'^ is the one that 
removes no corners {if X V' n), the next /i is the one that removes the highest 
corner, the next /i is the one that removes the second highest corner, etc. 



10.3 Proof of Theorem fnTD 

To provide motivation for the proof of Theorem 110.11 we begin by proving the 
recursive formula 

Theorem (Theorem [12]). For n > 3, |i?„| = 2n|i?„_i| - (n - l)2|i?„_2|- 

Proof. Viewing the elements of i?„ as rook matrices, i?„ consists of those el- 
ements having all O's in column n and row n (of which there are to- 
gether with, for each a G {l,...,n}, those having a 1 in position (a,n) (of 
which there are n|i?„_i| total), together with, for each a e {1, .... n — 1}, those 
having a 1 in position (n, a) (of which there are (n — l)|i?„_i| total). Counting 
the number of elements of i?„ in this way overcounts. For each {a, (3) with 
1 < Of, /3 < n — 1, every element with I's in positions (a, n) and (n, /?) (of which 
there are [n — l)^|i?„_2| total) gets counted twice. □ 

Now, let / e Ci?„ be given with respect to the {s} basis. We prove Theorem 



Theorem (Theorem 110. 1|) . For any representation p on Rn, if n > 3, we have 
the following factorization. 

n 

/(p)^5]p(T,) J2 fTMpi^)+Pi[n]) E f[nMP(') 

2 = 1 Sei?„_l sGRn-l 



71-1 

E 

z=l 



p(r 
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where [n] is the link (1)(2) . . . (n — /a(s) — f{As), and 

j.rpz ifn — le dom(s) 

[/(sT*) otherwise. 

Proof. Let n > 3. We have 3 types of elements s of 

• Type 1: Those for which s{n) = i for some 1 < i < n. 

• Type 2: Those for which both s(i) = n for some 1 < i < n — I and 
n ^ dom(s). 

• Type 3: Those for which both s{i) ^ n for all 1 < i < ?^ and n ^ dom(s). 

By the argument given in the proof of Theorem 12. 2i this counts all elements of 
Rn precisely once. 

Now, let 1 < j < n. View the "coset representative" Ti as a permutation 
matrix, and view the elements s € Rn as rook matrices. Multiplying any matrix 
X on the left by Ti simply moves row j of X to row j + 1 (for all j such that 
'i j ^ n — 1) and moves row n of X to row i. Thus, as s varies over Rn-i, TiS 
varies bijectively over {s S Rn ■ s{n) — i}. Therefore, we have 

n 

sGRn of Type 1 z^l s^Rn-i 

n 

= E E fTMpm)pis) 

n 

i=l seRn^i 

where /T.(s) = /(r.s). 

Similarly, multiplying any matrix X on the right by moves column j of 
X to column j + l{i<j<n~l) and moves column n of X to column i. Thus, 
as s varies over Rn-i, sT^ varies bijectively over {s e i?„ : s{i) — n}. So 

E fi')pi')- E fi^npi^n 

s^Rn:s(i)—n s^Rn-i 

= E fisnpis)pir). 

se-R„-i 

To ensure that we only count the elements of Type 2, we restrict our attention 
tol<i<n— 1, and we set the function values of the elements of Type 1 
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appearing in the above sum to 0: 



seRr^ of Typo 2 i=l s6i?,„_i 



= E 

i=l 



Pin, 



where 



if n — 1 € dom(s) (i.e. n € dom(sT*)) 

f{sT^) otherwise. 



Finally, 



E E /(N^)p(N^) 

sS-Rn of Typo 3 sei?.„_i 

Putting this all together, then, we find that for any representation p of _R„ and 
any n > 3, 

f{p)^ E fi^)pi^)+ E /(^)/'(^)+ E /(^)^(^) 

seRr^ of Type 1 of Type 2 sG-Rn of Type 3 

n 

^T.P^T,) E fTMp{^)+p{[n]) E /n(^)/^(^) 



i=l 

n-1 

+ E 

i=l 



s£R„ 



□ 



Concluding remarks 



The extension of FFTs to semigroups creates a new collection of interesting 
challenges. We remark that many of the ideas in this paper (such as those 
from [3T] and several of the results in this paper that follow from them) can be 
extended to general inverse semigroups. As with groups, FFT algorithms for 
inverse semigroups will vary from semigroup to semigroup, but a number of the 
underlying ideas are the same for any inverse semigroup FFT. These ideas can 
be cast in a general framework, which we intend to help guide the development 
of future FFTs, and this general framework is the subject of a paper currently 
in preparation. 
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Also, as mentioned previously, spectral analysis for the rook monoid involves 
projecting a function onto the isotypic subspaces of Ci?„, which can be accom- 
plished by means of an FFT, and examining the resulting projections. A variety 
of interesting issues arise in this analysis, and we have worked out a detailed ex- 
ample, consisting of partially ranked voting data on R^, to explain them. These 
results are the subject of a paper currently in preparation. 
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